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The extractions of nuclear charge radii from muonic atom spectroscopy for *°Ca and ?°°Pb are 
revisited to analyze the model dependencies induced by employing a Fermi-type charge distribution. 
For that, the charge densities, together with the corresponding muonic transition energies, calculated 
by the covariant density functional theory are used as a benchmark. The root-mean-square deviation 
of transition energies is calculated to quantitatively investigate the sensitivities of transition energies 
to the details of the two-parameter Fermi distribution. It is found that the second and fourth 
moments of the charge distribution can be extracted accurately from the muonic atom spectroscopy 
without much model dependencies, whereas the obtained two-parameter Fermi distributions cannot 
reproduce the details of the benchmarking charge densities and, in particular, its surface-diffuseness 
parameter cannot be determined accurately with the present experimental uncertainties on the 


muonic transition energies. 


I. INTRODUCTION 


Muonic atoms, where a negatively charged muon is 
captured by a nucleus, have been widely used to extract 
nuclear charge radii with high-precision for most stable 
isotopes [1-6], since the technique of X-ray spectroscopy 
from ji-mesonic atoms developed in 1953 [7]. In addition, 
nuclear electric quadrupole moments can also be accu- 
rately deduced from the observed hyperfine splittings of 
muonic X-ray transitions [8-10]. However, there exists a 
long-standing problem that limits the extractive accuracy 
and leads to the non-negligible theoretical uncertainty, 
when the experimental transition energies are analyzed, 
i.e., there exists a puzzling discrepancy between theoret- 
ical calculations and experimental data for fine-structure 
splittings [3, 11]. 

In the early analysis, the nuclear polarization (NP) 
correction is supposed to be the main source of such a 
discrepancy, due to the existence of challenging prob- 
lems in the calculation of the NP effect [3, 11]. Thus, 
considerable theoretical efforts have been made in re- 
cent decades to perform a rigorous calculation of NP cor- 
rection [12-16]. Employing state-of-the-art techniques, 
the most comprehensive calculations of NP corrections 
are performed recently by Valuev et al. [16] and it is 
concluded that the NP effect is unlikely responsible for 
the discrepancy. As a result, the other effects could 
play roles in explaining the discrepancy. On the other 
hand, note that the influence of model dependencies, in- 
duced by employing a Fermi-type charge distribution in 
the model-dependent analysis, is in general not exam- 
ined [3]. Although the influence of model dependencies 
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could be evaded in a model-independent combined analy- 
sis of muonic data and electron-scattering data as shown 
in Ref. [11], it is also important to investigate the contri- 
bution of model dependencies to the uncertainty of ex- 
tracted charge radii. For this purpose, a detailed analy- 
sis of model dependencies with respect to the Fermi-type 
charge distribution is performed in the present work by 
using the theoretical charge density distribution from the 
self-consistent mean-filed calculations as a benchmark. 

In recent decades, the covariant (relativistic) density 
functional theory (CDFT) has attracted extensive atten- 
tions on account of its successful description of many nu- 
clear phenomena [17-22]. For instance, it can well re- 
produce the isotopic shifts in the Pb [23] and Sn [24] re- 
gions, and naturally give the origins of the pseudospin 
symmetry in the nucleon spectrum [25] and the spin 
symmetry in the anti-nucleon spectrum [26, 27]. In ad- 
dition, it can provide a consistent description of cur- 
rents and time-odd fields, which play an important role 
in nuclear magnetic moments [28-32] and nuclear rota- 
tions [33, 34]. Based on the CDFT with the Bogoliubov 
transformation in the coordinate representation, the rel- 
ativistic continuum Hartree-Bogoliubov (RCHB) theory 
was developed [35, 36] to provide a proper treatment of 
pairing corrections and mean-filed potentials in the pres- 
ence of the continuum and has achieved great success 
in various aspects about exotic nuclei, such as provid- 
ing a microscopic and self-consistent description of halo 
in ''Li [35], predicting the giant halos phenomena in 
light and medium-heavy nuclei [37-39], and reproducing 
the interaction cross sections and charge-changing cross 
sections in sodium isotopes and other light exotic nu- 
clei [40, 41]. 

In this work, nuclear charge densities from the RCHB 
calculations are used as a benchmark to revisit the ex- 
tractive process of charge radii based on two-parameter 
Fermi (2pf) distribution for two doubly-magic nuclei, 


40Ca and ?°8Pb. In addition to the second moment, the 
fourth moment of nuclear charge distribution has recently 
attracted great attention for various applications, such as 
providing information on nuclear surface properties [42] 
and performing a successful study of the neutron-skin 
thickness [43]. In this paper, the sensitivities of muonic 
transition energies to the second and fourth moments are 
also investigated. 

The paper is organized as follows. In Sec. II, we in- 
troduce the RCHB method for solving the nuclear struc- 
ture and then construct the nuclear charge densities. In 
Sec. II A, we determine the parameters of 2pf distribu- 
tion by fitting the charge densities from RCHB calcu- 
lations in several ways and compare the corresponding 
transition energies with the data of RCHB. In Sec. IIIB, 
we determine the parameters by fitting the referenced 
transition energies, namely the results with the charge 
densities from RCHB calculations. In this way, we quali- 
tatively and quantitatively investigate the model depen- 
dencies and the sensitivities of transition energies to the 
second and fourth moments. Then, an analysis of the 
Barrett model is also made for further comparison. Fi- 
nally, summary and perspectives are presented in Sec. IV. 


II. THEORETICAL FRAMEWORK 
A. Nuclear charge density from RCHB theory 


In this study, the RCHB theory constructed with the 
contact interaction in the point-coupling representation 
between nucleons is adopted, where the conventional 
finite-range meson-exchange interaction is replaced by 
the corresponding local four-point interaction between 
nucleons. The details of the RCHB theory with point- 
coupling density functional can be found in Refs. [44, 45]. 
In the following, we introduce the theoretical framework 
briefly. 

Starting from the Lagrangian density, the energy den- 
sity functional of the nuclear system can be constructed 
under the mean-field and no-sea approximations. By 
minimizing the energy density functional with respect to 
the densities, one obtains the Dirac equation for nucle- 
ons within the relativistic mean-field framework [46]. The 
relativistic Hartree-Bogoliubov model provides a unified 
description of both the mean field and the pairing corre- 
lation. 

In the RCHB theory, the proton and neutron densities 
can be constructed by quasiparticle wave functions, 


pr(r) = X Vi (r)Ve(r) (1) 


ker 


with r € {p,n}. The nuclear charge density then includes 
the contributions from the point neutron density, the 
proton and neutron spin-orbit densities, and the single- 
proton and single-neutron charge densities, in addition 
to that from the point proton density [47-50]. The rela- 


tivistic nuclear charge density is finally written as 

Pe(r) = 3, [Per(r) + Wer(r)] , (2) 
where 
1 oo 
Pelr) = — f zpr (a) [gr(Ir — x|) — g(r + J] dz, (3) 
1 oo 
Werlr) == | 2We(2) [far(ln = al) = far(r +a)ldn. 
0 

(4) 
Here, p;(r) is the point nucleon density in Eq. (1) and 


W,(r) is the spin-orbit density given in Refs. [48, 50]. 
The functions g,(x) and for(x) are given by 


1/9, 
| Cora, (5) 
fer) =; fe! Far(q?) da, (6) 


in which Gg, and Fa, denote the electric Sachs and Pauli 
form factors of a nucleon, respectively. The following 
forms [48] are used in this study, 


1 
Grp’) = rera 

1 1 
Gen(g?) = az r2.q2/12)° (1+ r2 2/12)” 
Fæ) = Er 
ae ” 


with the proton charge radius r, = 0.8414 fm [51] and 
ri =r2,+ $(r2), where r2, = 0.9 fm? is the average of 
the squared radii for positive and negative charge distri- 
butions, while (r2) = —0.11 fm? [52] is the mean squared 
charge radius of a neutron. See reference [53] for more 


details. 


B. Dirac equation for muonic atom 


The Dirac equation for a muon reads 


[a-p+ 6M, + V(r) V(r) = ed(r), (8) 


where the eigenenergy e includes both the muonic atom 
energy levels E and reduced mass My, i.e., 


Mam, 


~ E+ M,, A ., 
ae Ma +m, 


with M, = (9) 


Note that Eq. (9) holds exactly only at the non- 
relativistic limit, while it is a good enough approxima- 
tion in the present calculations. The mass of muon 


My, = 206.7682830 m. is adopted [51] and nuclear mass 
Mz, is obtained from Ref. [54]. Since the electrostatic 
potential V(r) made by the atomic nucleus is spherically 
symmetric, the eigenvalue function can be written as 


_1(  iPu(r)¥fn(,9) 
Dnam(t) = — a (a €) Ym (0, ») 


where P„“(r) and Qn«(7) are its large and small compo- 
nents, respectively. 

The radial Dirac equation for one-muon system can be 
then written in the form of 


(of Ps) vail) (Ge) <2 (S28). 


r 


(10) 


where the mass term M, is subtracted on both sides of 
the equation. 
The electrostatic potential V(r) is obtained via [55] 


r yl oo 
V(r) = —4ra / p-(r’)— dr’ +f pe(r’)r’ ar! 
0 r r 
(12) 
where a is the fine structure constant. In the present 
work, the muonic energy levels are evaluated numerically. 


II. RESULTS AND DISCUSSION 


In order to investigate the influence of model depen- 
dency on muonic atom spectrum, nuclear charge densi- 
ties of Ca and ?°°Pb obtained from the RCHB calcu- 
lations are used as benchmarks. The relativistic density 
functional PC-PK1 [44], which provides one of the best 
density-functional descriptions for nuclear properties [56- 
63], is employed. The box size Ryox = 20 fm, the mesh 
size Ar = 0.1 fm, and the angular momentum cutoff 
Jmax = 19/2 ħ are used in the RCHB calculations. More 
numerical details can be found in Ref. [45]. 

The Dirac equation for the muon is solved using the 
generalized pseudospectral (GPS) method, whose pow- 
erful performance has been shown in Refs. [64-67]. The 
mapping parameters L = 0.01 and k = 3 are used in the 
GPS calculations. 


A. Fitting the nuclear charge density 


In general, the available experimental data of muonic 
atom spectrum are analyzed by means of a 2pf distribu- 
tion in the form of 


p 
pel") = Tepid Ins = e)/t]' = 


with the half-density-radius c and surface-diffuseness pa- 
rameters t, together with the normalization 


an | pe(r)r? dr = Z, (14) 
0 


where Z represents the proton number. The n-th mo- 
ment is given by 


4 co 
Rn = (r") = 5 | pel(r)r™*? dr. (15) 


In the following investigation, first of all, the parame- 
ters c and t of 2pf distribution are determined by fitting 
the RCHB charge densities in several feasible ways, with 
po determined by the normalization. Specifically, four 
2pf distributions are obtained here, i.e., 2pf(oı) is de- 
termined by minimizing the loss function o1, 2pf(Re2,01) 
and 2pf(R2,02) are determined by minimizing respec- 
tively gı and og with the second moment Rə constraint, 
and 2pf(R2,R4) is determined by constraining the sec- 
ond and fourth moments. The loss functions cı and oa 
are two alternatives that indicate the deviations in shape 
with the following forms: 


a= f (E-E) ar a) 
and 
o= | (er ran an 


As shown in Fig. 1, four 2pf distributions are compared 
with the charge density from the RCHB calculations. It 
can be seen from Fig. l(a) that the central densities of 
four 2pf distributions for 4°Ca are visibly different. The 
curve of 2pf(R2,R,) deviates furthest from the RCHB 
calculation, and the curve of 2pf( R2,02) is the next. In 
contrast, the curve of 2pf(o1) is closest in shape to the 
RCHB calculation. As for 2°°Pb in Fig. 1(b), it shows 
less differences between four 2pf distributions. Neverthe- 
less, the curve of 2pf(Ra,Rı) also shows the most devia- 
tion from the RCHB calculation. On the other hand, in 
the insets of Fig. 1 the logarithmic coordinates are taken 
to emphasize the difference of charge density at the sur- 
face. One can see from the insets that the distributions 
of 2pf(Ra,Rı) are the closest ones to the RCHB charge 
density distributions among four distributions at the nu- 
clear surface, i.e., r > 6 fm for Ca and r > 8 fm for 
208Ph, In contrast, the charge distributions 2pf(cı) are in 
the worst agreement with the RCHB calculations. With 
regard to the rest of two 2pf distributions, the curve of 
2pf(R2,02) is closer to (further away from) the curve of 
RCHB than 2pf(R2,01) for ?%Ca (?°8Pb). 

The fitted parameters c and t of four 2pf distributions, 
together with the corresponding second and fourth mo- 
ments and charge radii r., as well as two loss functions 
cı and oa are shown in Table I. Since the results of 
2pf(oı) are obtained without the second-moment con- 
straint, there exists a visible difference in the second mo- 
ment between 2pf(c1) and the others, i.e., around 0.7 fm? 
for °Ca and 0.2 fm? for 2°°Pb. Moreover, it can be seen 
that the parameters determined in 2pf(o1), 2pf(Ra,cı), 
and 2pf( R2,02), without the fourth-moment constraint, 
generally overestimate the fourth moments. 


= 2pf (R; 62) i 
---- 2pf (R, Ry) 1 


r (fm) 


r (fm) 


FIG. 1. Comparison of charge densities from the RCHB theory with those of 2pf distributions for (a) “°Ca and (b) 7°°Pb. The 


2pf parameters are determined by fitting the RCHB charge densities in four different ways, i.e., 2pf(o1 


and 2pf(R2,R4). See texts for details. 


), 2pf(R2,01), 2pf(R2,02), 


TABLE I. The parameters c and t, the second and fourth moments Ra and R4, and the charge radii re, as well as the values of 


two loss functions c1 and os, for *°Ca and ?®Pb in four 2pf models. Note that the numbers that exactly reproduce the RCHB 
values are highlighted in bold form. 
Ca EP 
2pf(oı) 2pf( R2,01) 2pf( R2,02) 2pf(R2,R4) 2pf(o1) 2pf( R2,01) 2pf( R2,02) 2pf(R2,R4) 
c (fm) 3.43431 3.51641 3.63210 3.72839 | 6.65823 6.66207 6.64228 6.70918 
t (fm) 2.82022 2.55115 2.41134 2.28476 | 2.33968 2.26797 2.31614 2.14837 
Ra (fm?) 12.7666 12.0755 12.0755 12.0755 | 30.5160 30.3103 30.3103 30.3103 
Ra (fm*) 256.054 220.280 214.499 209.699 |1195.877 1174.858 1178.577 1166.040 
re (fm) 3.5730 3.4750 3.4750 3.4750 5.5241 5.5055 5.5055 5.5055 
1 (10~°fm~*)|32.1726 44.8123 90.7370 195.203 | 9.71232 10.42702 11.3099 15.5466 
2 (10~°fm~*)}19.1357 18.4373 9.36914 16.1918 | 10.4874 13.9824 12.6133 29.1877 


TABLE II. Muonic transition energies (keV) calculated by using RCHB charge density, and the differences between them and 
the results of four 2pf distributions, for *’Ca and 7°°Pb. Numbers in parentheses represent the power of ten. 


Transition ErRcus Ca, AEspt — AERcHB ERcuB 208Pp, ABsapt — AERcHB 

(keV) (4°Ca)|2pf(o1)  2pf(Ra,0ı) 2pf(R2,02) 2pf(Re,R1)| (?°8Pb)|2pf(o1)  2pf(R2,01) 2pf(Re,o2) 2pf(Re,Ra 
2p3/2-181/2 | 778.229 | —2.58 1.82(-1) 8.47(-2) 3.69(-3)|5917.293] —8.77 4.84 7.06 —4.51(—1 
2p1/2-181/2 | 776.766 | —2.58 1.82(—1) 8.47(—2) 3.68(—3)]5735.444|—8.24 4.67 6.82 —4.57(—1 
3d3/2-2p1/2| 157.439 | —8.03(—3) —1.06(—3) —4.46(—4) 7.14(—5)|2619.561 | —2.69 7.98(—4) —9.02(—2) 2.15(—1 
3d5/2-2p3/2 | 156.124 | —5.43(—3) —1.11(—3) —4.74(—4) 5.95(—5)|2479.793 | —2.20 1.78(—1) —3.40(—1) 2.12(—1 
3d3/2-2p3/2 | 155.976 | —5.43(—3) —1.11(-3) —4.74(-4) 5.96(—5)|2437.712| —2.16 1.72(—1) —3.31(—1) 2.09(—1 
5 fz/2-3d5/2| 79.851|—1.96(—6) —6.14(—7) —3.05(—7) —5.66(—8) |1358.683|—5.67(—2) —1.88(—2) —2.75(—2) 1.87(—3 
5 fs/2-3d5/2| 79.835|—1.95(—6) —6.05(—7) —2.96(—7) —4.78(—8) |1354.023|—5.64(—2) —1.87(—2) —2.74(—2) 1.87(—3 
4 fz/2-3d5/2| 54.599|—1.96(—6) —6.07(—7) —2.99(—7) —5.04(—8)| 932.086|—5.66(—2) —1.88(—2) —2.74(—2) 1.87(—3 
4 fs/2-3d5/2| 54.568 |—1.95(—6) —6.06(—7) —2.98(—7) —4.95(—8)| 923.018|—5.63(—2) —1.87(—2) —2.73(—2) 1.87(—3 
4ds/2-3p1/2| 55.074|—2.81(—3) —3.70(—4) —1.55(—4) 2.53(—5)| 914.494|—7.84(—1) 4.73(—2) 4.00(—2) 6.57(-2 
4ds/2-3p3/2| 54.705|—1.91(—3) —3.90(—4) —1.66(—4) 2.13(—5)| 885.170|—6.94(—1) —2.38(—2) —6.41(—2) 7.31(—2 
4d3/2-3p3/2| 54.643 |—1.91(—3) —3.90(—4) —1.66(—4) 2.13(—5)| 867.687|—6.75(—1) —2.12(—2) —5.96(—2) 7.15(—2 


In Table II, the muonic transition energies calcu- 
lated by using the RCHB charge density, and the dif- 
ferences between them and the results calculated by us- 


ing the 2pf distributions, for *°Ca and ?°°Pb, are shown. 
Twelve muonic transition energies for *°Ca and ?°8Pb 


are given, namely, 2p3/2-151/2, 2pı/2-151/2, 3d3/2-2p1/2, 


3d5/2-2p3/2, 3d3/2-2p3/2, 5 f7/2-3ds/2, 5f5/2-3ds5/2, 4f7/2- 
3d5/2, 4f5/2-3d5/2, 4d5/2-3p1/2, 4d5/2-3p3/2, and 4d3/2- 
3p3/2, involving fourteen levels. 

In comparison, the results of 2pf(o1) are in the worst 
agreement with RCHB calculations and the results of 
2pf(R2,R4) are in the best agreement with RCHB calcu- 
lations. For example, the 2pf(a,) gives 2.58 keV deriva- 
tion from RCHB for 2p3/9-1s1/2 in *°Ca. After con- 
sidering the second moment constraint Ra with differ- 
ent loss function minimized, 2pf(R2,01) and 2pf(Re2,02) 
give the derivations —0.182 and —0.0847 keV, respec- 
tively. Furthermore, with both second- and fourth- 
moment constraints, the derivation from RCHB calcu- 
lations decreases to —3.65 x 1073 keV. A similar conclu- 
sion for other transition energies in Ca and ?08Pb is 
also observed in Table II. In view of the comparison of 
charge distribution between RCHB calculations and four 
2pf distributions shown in Fig. 1, it is concluded that 
the surface of nuclear charge density plays a much more 
important role in the muonic atom spectrum than the 
center. Furthermore, it implies the considerable sensitiv- 
ities of the transition energies to the second and fourth 
moments. A quantitative discussion is given in the next 
section. 


B. Fitting the muonic transition energies 


In this section, the procedure to extract the charge 
radii from the muonic transition energies with the 2pf 
distribution will be discussed. For the sake of discus- 
sion, the root-mean-square deviation (RMSD) for muonic 
transition energies is defined as 


1/2 
= È D (ar? = az)" , (18) 


[7 


where E? Pf and ERCHB are the i-th transition energy cal- 
culated with the 2pf charge distributions and the RCHB 
charge distributions, respectively, and N is the total 
number of transitions considered. For convenience, four 
transitions with the largest transition energies, i.e., 2p3/2- 
1s1/2, 2p1/2-181/2, 3d3/2-2p1/2, and 3d5/2-2p3/2, are con- 
sidered in Eq. (18). 

In panels (a) and (d) of Fig. 2, the trend of 6 with 
the parameter t changing from 1 to 3.5 fm is shown, 
where the parameter c is determined by minimizing the 
RMSD ô for each t. A minimum point of ô can be 
seen in the panels. Here, we use the shorthand writ- 
ing 2pf(AE) to denote the best-fit 2pf distribution which 
yield a minimum value of 6. The minimum values of 
ô read 6 > 2.0 x 1078 keV with t = 2.2950 fm and 
c = 3.7214 fm for Ca and 6 > 0.0031 keV with 
t = 2.1821 fm and c = 6.6977 fm for ?0®Pb. In addi- 
tion, ö versus the second moment (r?) of the 2pf dis- 
tribution with the same parameters as Figs. 2(a) and 
2(d) is also displayed in Figs. 2(b) and 2(e). The sec- 
ond and fourth moments of the best fit are determined 


as the point of minimum 6, namely (r°) = 12.0778 fm?, 
(r*) = 210.145 fm* for Ca, and (r?) = 30.3225 fm?, 
(r*) = 1169.38 fm? for SPb. Meanwhile, the referenced 
second moments from the RCHB calculations are also 
shown in the panels. The relative deviations in charge 
radii read |v 12.0778 — v12.0755|/v 12.0755 = 0.01% for 
40Ca and 0.02% for ?°°Pb. Note that the source of the 
deviation here is the model dependency. From the above 
numbers, it can be seen that the influence of the model 
dependency on the extracted charge radii is rather small. 

In practice, the accuracy of the extracted charge radii 
is limited by the total uncertainty of transition energies 
from experimental measurements, together with various 
additional corrections, such as vacuum polarization, rel- 
ativistic recoil, self-energy, and nuclear polarization cor- 
rections [68-71]. Therefore, the determination of uncer- 
tainty of the extracted charge radii is in general a chal- 
lenging task. In the present theory-to-theory benchmark- 
ing study, the dependence of the uncertainty in the ex- 
tracted charge radii Ar. on the uncertainty in the tran- 
sition energies can be quantitatively investigated. The 
RMSD ô in the present calculations can be regarded as 
the total uncertainty of transition energies, and the dif- 
ference of charge radii Ar, = r?Pf — rRCHB can be used 
to indicate the uncertainty of the extracted charge radii. 
For this purpose, the variations of 6 versus Ar, for Ca 
and ?°8Pb are shown in Figs. 2(c) and 2(f), respectively. 
The points at Ar, = +0.01 fm are marked and the cor- 
responding values of 6, i.e., 6 © 0.0017 keV for Ca and 
ô = 1 keV for ?0®Pb, are illustrated. In fact, it is quite 
challenging to make the total uncertainty of transition 
energies less than 0.0017 keV for Ca. However, the un- 
certainty less than 1 keV for ?°8Pb is available [3]. As 
a result, the uncertainty of the extracted charge radii 
less than 0.01 fm for ?°°Pb can be obtained. It can be 
deduced that the charge radii for the heavy nuclei can 
be extracted from muonic atom spectroscopy more accu- 
rately than for the light nuclei. 

In order to study the variation of 6 versus the fourth 
moment, the influence of the second moment should be 
avoided. In Fig. 3, the steep tendency of RMSD 6 with 
the fourth moment is illustrated, where parameter t is 
changed from 1 to 3.5 fm, and c is determined accord- 
ing to the additional constraint of the second moment. 
The fourth moment of the RCHB charge distributions 
is marked and a small deviation from the corresponding 
value at the point of minimum 6 can be seen. The rel- 
ative deviations of the fourth moment are only around 
0.2% and 0.3% for Ca and ?08Pb, respectively. There- 
fore, it can be considered that the prediction of the fourth 
moment from the muonic atom spectrum is reliable. 

In addition, the analysis using the Barrett model [1] 
is also performed. The Barrett equivalent radius Rķa is 
defined by 


3 Rra 
af pe Pr dr = irre). (19) 
ka J0 


It can be iteratively calculated from the charge density 
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FIG. 2. The variation of ô for *°Ca and ?°8Pb: (a) and (d) 6 versus parameter t of 2pf distribution, where parameter c are 
determined to minimize 6; (b) and (e) 6 versus the second moment of 2pf distribution, where the black dashed lines indicate the 
second moment of RCHB charge distributions; (c) and (f) 6 versus the difference of charge radii between the 2pf distributions 


and RCHB calculations, i.e., Are = r?Pf — rRCHB 


without resolving the Dirac equation. The critical sensi- 
tivity of transition energies in muonic atoms to the Bar- 
rett equivalent radius has been illustrated [72, 73]. Thus, 
the Barrett model provides a particular perspective to 
compare these 2pf distributions with the fitted param- 
eters to the referenced charge density from the RCHB 
calculations. The comparisons of the Barrett equivalent 
radius for four 2pf distributions and the RCHB calcula- 
tions versus k for Ca and ?P®Pb are shown in Fig. 4. 
The values of a are in general determined by the fitting 
of the differences of the potential that muon in state a 
makes and the potential that muon in state b makes, i.e., 
fav(r) = Va(r) — Vier) [3, 72]. The muon potential in the 
state a is defined by [1] 


nf Pal Higi EL igh (20) 


where P,(r) and Qa(r) are the large and small compo- 
nents in Eq. (11), respectively. In the present calculations 
of the Barrett equivalent radius, the values of a are taken 
as 0.065 fm! for "Ca [2] and 0.1415 fm! for ?°8Pb [3]. 
Different from the practical analysis in experimental ref- 


, where the black dashed lines represent Ar. = +0.01 fm. 


erences, e.g., [3], the results of benchmark, i.e., “RCHB”, 
are displayed in the form of smooth curves instead of sev- 
eral dots that correspond to transition energies, because 
they are evaluated from the RCHB charge density instead 
of the transition energies. It can be seen that the results 
of 2pf(R2,R4) are in excellent agreement with the results 
of RCHB calculations. In contrast, the results of 2pf(o1 ) 
considerably deviate from the results of RCHB calcula- 
tions. As for the rest of two 2pf distributions, Fig. 4 
shows that the result of 2pf( R2,02) is in a better (worse) 
agreement with the result of RCHB than 2pf(Ra,cı) for 
40Ca (208Pp), which are consistent with the comparison 
of transition energies displayed in Table II. 

Furthermore, the results of 2pf(AE), i.e., the best-fit 
2pf distributions for minimizing the RMSD ö mentioned 
above, are also compared in the Fig. 4. It can be seen 
that these results are also in excellent agreements with 
the results of RCHB calculations. It ensures that the 
parameters fitted to the transition energies give a satis- 
factory agreement. 

In order to perform a more practical investigation, we 
will then focus on the weighted RMSD, which is defined 
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FIG. 3. The root-mean-square deviation 6 as a function of 
the fourth moment R4 of 2pf distribution for (a) Ca and 
(b) SPb. Different from Fig. 2, an additional constraint 
is imposed to keep the second moment constant, i.e., R2 = 
12.07784 fm? for “°Ca and 30.32254 fm? for ?°°Pb. The black 
dashed lines indicate the data from RCHB calculations. 
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where n; represents the relative deviations for i-th tran- 
sition energy. According to one of the latest experiments 
for muonic 1°*Pd [74], the relative uncertainties of experi- 
mental transition energies are approximately 6.47 x 1075, 
9.81 x 1075, 3.47 x 1074, and 1.20 x 1074 for transitions 
2p3/2-181/2, 2P1/2-181/2, 3d3/2-21/2, and 3d5 /2-2p3/2, re- 
spectively. Thus, the corresponding n; in Eq. (21) are 
taken as 1074, 1074, 3x 1074, and 3 x 1074, respectively. 

In Fig. 5, the confidence region that satisfies the condi- 
tion x < 1 for ?P®Pb is displayed. Parameter t is changed 
from 1 to 3.5 fm with a step 0.01 fm, and c is determined 
to minimize x for each t. Then two critical points could 
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FIG. 4. The Barrett equivalent radius Rpa for (a) Ca and 
(b) SPb. The values of a in Eq. (19) are taken as 0.065 fm! 
for “°Ca [2] and 0.1415 fm! for *°8Pb [3]. 


be found and the confidence region 2.05 < t < 2.31 is 
obtained in Fig. 5(a). This implies that the relative un- 
certainty of t is around (2.31 — 2.05)/(2.31 + 2.05) = 6%, 
while the relative uncertainty of c reaches around 0.6%. 

In Fig. 5(b), the Barrett equivalent radius of the 2pf 
distributions with parameters corresponding to two crit- 
ical points are compared with the results of RCHB cal- 
culations. It is of interest that three curves intersect ap- 
proximately at a point around k ~ 2.25, and the curve 
of RCHB is sandwiched between two 2pf curves with k 
in the range from 1 to 4.5. 

Furthermore, in Fig. 6, the elliptic region of 2pf pa- 
rameters which satisfy the condition x < 1 is shown. The 
corresponding charge radii are in the range from 5.5004 
to 5.5129 fm. Thus, the uncertainty of the extracted 
charge radii is only around 0.005 fm, and the relative un- 
certainty is around 0.11%. As a result, it is considered 
that the extraction of charge radii from the muonic atom 
spectrum for ?°8Pb is reliable. 

Meanwhile, there are several interesting observations 
from the figure. First of all, it shows the contour of the 
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FIG. 5. Determination of parameters in 2pf model for ?°°Pb. 
(a) The weight RMSD x versus t, where c is searched to min- 
imize x. Two critical points satisfying condition x < 1 are 
indicated. (b) Comparison of the Barrett equivalent radius 
between the RCHB calculations and the 2pf distributions with 
respect to two critical point. 


rms radius in the c-t plane. One can see that the spread 
of the elliptic region is almost parallel to the contour of 
the rms radius. Combined with the comparison of rela- 
tive uncertainty between the parameters c and t and the 
extracted charge radius as previously discussed, it is con- 
firmed that the transition energy is more sensitive to the 
rms radius than to the individual parameter c or t. See 
also the discussions in Ref. [74]. Second, the contour of 
the referenced rms radius, i.e., ae is shown in Fig. 6. 
In addition, in the elliptic region, the point with the min- 
imum y and the line with the c determined to minimize 
x for each t is also shown. They correspond to the red 
star symbol and the blue curve in Fig. 5(a), respectively. 
Keeping in mind that, in this figure, the green dashed 
line (representing the determined c values by minimiz- 
ing x for each given t) would be identical with the red 
dot-dashed line (representing the contour of the bench- 
marking charge radius), if there is no model dependency 
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FIG. 6. Determination of parameters in 2pf model for ?°®Pb. 
The elliptic region of 2pf parameters c and t for x < 1 is 
denoted with the blue shadow area, where the red star symbol 
shows the point with the minimum x and the green dashed 
line shows the determined c values by minimizing x for each 
given t. The contours of the rms charge radii (0.01 fm interval) 
are illustrated with the black solid lines. Additionally, the red 
dot-dashed line represents the contour of the benchmarking 
charge radius. 


on extracting the charge rms radius. In contrast, it is 
seen from the figure that the slopes of these two lines are 
slightly different from each other. Precisely, the slopes 
are —0.32 and —0.38 for the former and latter, respec- 
tively. Such a difference can be used as an indicator of 
the model dependency. 

Before ending this section, it is also important to men- 
tion that one cannot constrain the charge radius of 4°Ca 
at all, by taking the same relative deviations n;. This 
is because, to extract the charge radius at the level of 
Ar. = +0.01 fm, the required accuracy of the transi- 
tion energies is much higher in the case of Ca, as we 
compared quantitatively in Figs. 2(c) and 2(f). 


IV. SUMMARY AND PERSPECTIVES 


In this work, the influence of model dependencies on 
the extraction of charge radii from the muonic atom spec- 
trum and the sensitivities of muonic transition energies 
to the second and fourth moments of charge distribu- 
tion are investigated by taking 4°Ca and ?°8Pb as exam- 
ples. The charge densities from the RCHB calculations 
are used as a benchmark. The parameters c and ¢ in 
the 2pf model are determined in several ways to fit the 
referenced charge densities, and the corresponding tran- 
sition energies are compared with the data from RCHB 
calculations. The sensitivities of transition energies to 
the second and fourth moments are illustrated. 

In detail, first of all, the variations of RMSD 6 of 
muonic transition energies versus parameters c and t are 


investigated. By changing t from 1 to 3.5 fm and deter- 
mining c iteratively to minimize 6 for each t, the steep 
trend of ö versus the second moment is clearly shown. 
Meanwhile, the variations of 6 versus Ar, are also shown 
to investigate the dependence of the uncertainty in the 
extracted charge radii on the uncertainty in the transi- 
tion energies. On the one hand, it is found that the to- 
tal uncertainty of transition energies should be less than 
0.0017 keV for *%Ca and 1 keV for 2°°Pb to make the 
uncertainty of extracted charge radii less than 0.01 fm. 
This illustrates that the charge radii of heavy nuclei could 
be extracted more accurately from muonic atom spec- 
troscopy than those of light nuclei. On the other hand, 
when the charge radii are extracted by minimizing 6, it 
is found that the relative deviations of extracted charge 
radii are around 0.01% for 4°Ca and 0.02% for ?08Pb, 
compared with the referenced charge radii. Therefore, it 
is considered that the model dependency only slightly in- 
fluences the extraction of charge radii. Second, the vari- 
ation of 6 versus the fourth moment of 2pf distribution 
is displayed with the constraint of the second moment. 
The sensitivities of transition energies to the fourth mo- 
ment can be clearly seen. Finally, the analysis with the 
Barrett model is performed. The Barrett equivalent radii 
obtained by various 2pf distributions are compared with 
the referenced results. It is clearly seen that the 2pf dis- 
tribution determined by the transition energies 2pf(AE) 
and that with the second- and fourth-moment constraints 
2pf( R2, R4) reproduce the benchmarking Barrett equiv- 
alent radii excellently, whereas those without the fourth- 
moment constraint 2pf(Ra, 01) and 2pf(Ra, o2) provide 
less satisfactory results, and that only focusing on the 
density distribution alone 2pf(oı) shows the result con- 
siderably deviating from the benchmark. 

As a step further, the analysis with respect to the 
weighted RMSD x with the relative deviations n; taken 
from the recent experimental data is performed, by tak- 
ing ?0®Pb as an example. The confidence region for pa- 
rameters c and t, satisfying the condition x < 1, is ex- 
hibited in the c-t plane. The small uncertainty of charge 
radius, i.e., 0.005 fm, emphasizes the reliability of the 
extracted charge radii from the muonic atom spectrum. 
Meanwhile, it is considered that the difference of slopes 
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of the two lines—one represents the determined c values 
by minimizing y for each given t and the other repre- 
sents the contour of the benchmarking charge radius— 
in the c-t plane can be used as an indicator of the 
model dependency. However, it is also found that the 
obtained 2pf distributions cannot reproduce the details 
of the benchmarking charge densities. In particular, its 
surface-diffuseness parameter t cannot be determined ac- 
curately with the present experimental uncertainties on 
the muonic transition energies. 


In the present theory-to-theory benchmarking study, 
we have analyzed the model dependencies and deter- 
mined the uncertainty of extracted charge radii from 
muonic atom spectroscopy based on the 2pf distribu- 
tion. It is also expected to see how much information on 
the nuclear structure can be extracted from the muonic 
atom spectroscopy using a model-independent distribu- 
tion, such as the Fourier-Bessel series expansion [75]. 
This question would be of great interest in our future 
studies. 
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